Diffusion limited friendship network: A model for six degrees of separation 
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A dynamic model of a society is studied where each person is an uncorrelated and non-interacting 
random walker. A dynamical random graph represents the acquaintance network of the society 
whose nodes are the individuals and links are the pairs of mutual friendships. This network exhibits 
a novel percolation like phase transition in all dimensions. On introducing simultaneous death and 
birth rates in the population we show that the friendship network shows the six degrees of separation 
for ever after where the precise value of the network diameter depends on the death/birth rate. A 
SIS type model of disease spreading shows that this society remains always healthy if the population 
density is less than certain threshold value. 

PACS numbers: 05.10.-a,05.40.-a,05.50.+q,89.75.Hc 
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Though at present the human population of the world 
has attained a very large size, more than 6 x 10 9 precisely, 
it is believed that any two randomly selected persons in 
the world are connected by a short chain of intermediate 
acquaintances, typically of length 6. This phenomenon 
is referred as the "six degrees of separation". The idea 
originated from the famous letter distribution experiment 
of Milgram in 1960s lj. Since then any network of N 
nodes is said to display six degrees of separation if its 
diameter is small and grows at most as logiV 0,0- 

Most human communications, especially the informa- 
tion exchanges, take place directly between individu- 
als when they are at close proximity to one another. 
The spread of news, rumors, jokes and fashions all take 
place by communications among individuals. More im- 
portantly the infectious diseases also spread by person- 
to-person contact and the structure of network of such 
contacts has important effects on the nature of the epi- 
demics. Naturally the speed of spreading in general is 
faster for a network with small diameter. 

There are important models of the social networks like 
Small-world network (SWN) that displays the six degrees 
of separation 0. Also the process of spreading of epi- 
demics is modeled by a susceptible-infected-susceptible 
(SIS) model in which a non-equilibrium phase transi- 
tion takes place from a healthy society to an infected 
society at a critical value of the infection probability 

All these models of the social networks as well as for 
the spreading of the infectious diseases consider a static 
picture of the society. More precisely static individu- 
als are positioned at the nodes of certain graphs and a 
person interacts with only a fixed set of neighbours de- 
termined by the degree of the node. Where as in actual 
society the number of acquaintances of a person increases 
with time. Everyday a person goes to office, market, the- 
aters, clubs etc. and therefore gets acquainted with other 
people who were unknown to him. By the same move- 
ments a person becomes exposed to infections by others 
or transfers his own infection to others. Again not all 
the friends of an infected person has the chance of get- 



ting infection, only those friends who come close to this 
person has the risk of infection. In this paper we study 
this basic property of a dynamical society where indi- 
viduals are not static objects but move continuously and 
therefore comes in contact with other people. To make a 
simple model we have considered the diffusive motion of 
the individuals and modeled the society by a set of ran- 
dom walkers. Specifically in our model (i) unlike static 
models the number of acquaintances of a person evolves 
with time (ii) irrespective of how many friends an in- 
fected person has, he may infect only those friends who 
come to his close proximity, this is unlike to the ordi- 
nary SIS type models. What we achieve are: (i) with 
the introduction of a death/birth rate the society indeed 
shows the six degrees of separation effect (ii) there is a 
threshold density of population, below which the society 
is always healthy (iii) a very interesting theoretical obser- 
vation that the associated dynamical random graph has 
a non-trivial dimension dependent critical behaviour. 

Over last few years it is becoming increasingly evident 
that highly complex structures of many social |g, bio- 
logical electronic communication systems jlUE^ 
etc. can be modeled by simple graphs. Erdos and Renyi 
studied the well known random graphs (RG) of N nodes 
where each pair of nodes is connected with a probabil- 
ity p and the grap h shows a continuous phase transition 
at p c = 1/N [13. Scale- free networks (SFN) are char- 
acterized by the power law decay of the nodal degree 
distribution function: -P(fc) ~ fc~ 7 . Two very impor- 
tant networks in electronic communication system like 
World Wide Web 01 an d the Internet are observed 
to possess the scale-free property. Barabasi and Albert 
(BA) proposed a model for a growing SFN where nodes 
are linked with the preferential attachment probability 
3 E|- Other routes, e.g., static and quasistatic 
17] models to obtain SFNs are also studied. Assigning 
a Hamiltonian correlations are studied in the optimized 
networks keeping biological networks in mind |l8| • S WNs 
with random walkers capable of making long distance 
jumps are studied in [lflj . 

In our model each member of the society executes a 
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FIG. 1: Plot of the scaled characteristic times T c (opaque 
symbols) and T' c (filled symbols) vs. the system size L for 
three different populations: TV = 32 (circle), 64 (square) and 
128 (triangle). 



simple uncorrelated and non-interacting random walk on 
a regular lattice. Initially the population of TV persons 
are released on the square lattice of size L x L at randomly 
selected positions. The system then starts evolving with 
time. At each time step each person makes a jump to 
one of its neighbouring lattice positions with equal prob- 
abilities. Each person represents a node of the growing 
acquaintance network and a link is established between 
two nodes the moment the corresponding pair of persons 
come in contact to each other at the same position and 
at the same time [20j. Gradually the number of links 
among the individuals grow. Thus the set of TV nodes 
and the set of links among these nodes define our net- 
work called as the Diffusion limited friendship network 
(DLFN) where as the associated graph is referred as the 
Dynamical random graph (DRG). 

All persons which are at the same lattice site immedi- 
ately become friends and the associated subgraph with 
these people become a clique. At each time many such 
cliques are formed at different sites. All these cliques re- 
main for ever, they never get destroyed, moreover they 
grow in sizes as time proceeds. At the early times, the 
number of links is small and the DRG has many differ- 
ent isolated components of different sizes. The size of a 
component is determined by its number of nodes and the 
giant component has the largest size. The giant com- 
ponent not only grows by including new nodes into it 
but also by the process of merging equally large compo- 
nents. After some slow initial growth the giant compo- 
nent grows very fast and its size become proportional to 
TV. This behaviour is just like the threshold phenomenon 
in a continuous phase transition e.g., what happens in a 



FIG. 2: The scaling of the order parameter ip(p,N) at the 
critical region p a (N) in the DRG: N=64 (circle), 128 (square) 
and 256 (triangle). 



random graph |13|. The whole DRG ultimately reaches 
the limiting stage of a giant TV-clique when each node is 
linked to all other nodes. 

We first characterize DRG to compare with RG. DRG 
has two characteristic time scales. T c measures the time 
required for the phase transition and is observed to vary 
like L z /N a . In a mean- field limit when the density 
p = N/L d is small this variation is estimated in the fol- 
lowing way: If a person randomly walks a linear distance 
R in (i-dimension in time T c then R ~ and therefore 
around L d / R d such <i-dimensional spheres are needed to 
cover the volume of size L which is TV itself. This gives 
% ~ L^/TV 1 /^ i.e., z=l/n and a = l/fid in general 
and therefore z — 2 and a = 1 for ordinary random walks 
(/j = 1/2) in two dimension. Fig. 1 shows the scaled plot 
of % for different L and TV values and a good collapse of 
the data is observed but for a w 0.89 and z w 2.25. We 
believe the difference in the exponents from the mean- 
field values are due to finite size of the system. At a 
second characteristic time T' c ~ L z /N a the DRG be- 
comes an TV-clique where z 1 and a' are estimated as 2.22 
and -0.33 respectively for /i = 1/2 and d=2. The posi- 
tive value of a and the negative value of a' are consistent 
with intuition: for a fixed L but with increasing TV, less 
number of steps per person are necessary for the giant 
component to include all nodes, but a larger number of 
steps are required to form the TV-clique. The values of z 
and z 1 are likely to be the same. 

The link density p(t, TV) at a time t for an TV node 
network is defined as the ratio of the number of links to 
its maximum possible number TV(TV — l)/2. Numerically 
we find the following scaling form: 

P (t,L)~T(t/L z ) (1) 
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FIG. 3: The fluctuating diameters V(N,t) and V(N,t) of 
the DLFN (g=0.43, iV=64, L=64) plotted in the upper and 
the lower curves with time. The data is averaged over 100 
configurations and have steady averages of 6.00 ± 0.05 and 
2.87 ±0.05. 



FIG. 4: The variation of the activity A(r) of the SIS model 
on DLFN with the infection probability r for p = 1.2 and it 
vanishes at r c . The inset shows the variation of the critical 
infection probability r c (p) with population density p. 



where the scaling function T{x) ~ x a and a and z are 
approximately found to be 0.89 and 2.25 again. 

The order parameter ij)(p,N) of this transition is the 
average fraction of nodes in the giant component for a 
link density p. The critical link density at the transition 
point p c is defined by tp(p C7 N) = 1/2 and is observed to 
vary with N as: p c = b/N a with b « 1.28 and a w 0.89 as 
before. As the mean-field calculation gave a = 1/fJ.d, we 
see that only the ordinary random walks in two dimen- 
sions with p c — l/N correspond to the random graphs 
[l3j | but for other walks with different [i and in differ- 
ent dimensions p c (N) have non-trivial dimension depen- 
dence. A scaling plot for the order parameter is shown 
in Fig. 2 where we plot 4>{p, N)N-P/ V vs. \p-p c ]N 1 / v . 
An excellent collapse of the data shows that the order 
parameter has the following scaling form: 

i;(p,N)^N^g[(p- Pc )N 1 / 1 '} (2) 

with v k 0.8 and « 0.02. 

The topological distance between a pair of nodes is the 
number of links on the shortest path connecting them and 
the diameter is the maximum of such paths. The aver- 
age diameter T>{N) is measured over many independent 
configurations. The configuration average of the mean 
distance between an arbitrary pair of nodes is denoted 
by T)'(N). As the system evolves, both measures first in- 
crease with time, reach their maxima and then decrease 
very slowly, finally saturates to a fixed value for a long 
time. The maximum of the diameters occur at the char- 
acteristic times T c . As expected the nodal degree distri- 
bution of the giant component at the transition point is 
a Poisson distribution similar to RG, since there is no 



preferential link attachment probability in this model as 
in the scale-free networks. 

The network described so far has a major drawback 
that it assumes each individual as immortal. As a con- 
sequence the DRG becomes an iV-clique at time T' c . To 
make our model more realistic, we therefore introduce 
a probability of death and birth in the population but 
with equal rates to keep the population conserved. More 
precisely at each time step only one randomly selected in- 
dividual is killed with a probability q. As a consequence 
all links associated with the node representing this indi- 
vidual are immediately deleted. This may result the frag- 
mentation of the particular component of the dynamical 
graph which belonged this node. A fresh determination 
of the different components of the DLFN especially the 
giant component is done immediately before the system 
proceeds to the next time step. At the same time we as- 
sume that a fresh individual has taken birth at the same 
position of the dead individual so that the population 
conservation is maintained. 

When an individual dies, the deletion of all his links 
may severely affect the distribution of distances between 
all pairs of nodes in the system. In fact it is expected 
that in general the distance between an arbitrary pair 
of nodes should increase due to the death of an indi- 
vidual, which thus enhances the values of T>{N,t) and 
T>'(N,t). On the other hand the newly born individ- 
ual also start diffusing in the system and start building 
up links of acquaintances with other individuals of the 
network. Therefore the magnitudes of the diameters de- 
crease again. As a result the net effect of the two com- 
petitive processes of simultaneous death and birth of in- 
dividuals is to make the diameters fluctuate around their 
steady averages whose magnitudes must depend on the 
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rate q of death/birth processes. In Fig. 3 we show the 
time variation of these diameters and for L — 64, TV = 64 
the diameter T)(q, N, L) has a value very near to 6 for q = 

0. 43 where as the T>'(q, N, L) is around 2.85. As expected 
the diameters increase with decreasing q. 

Finally we study a susceptible-infected-susceptiblc 
model on the DLFN. At any time a lattice site may be 
occupied by a number of persons. If at least one of them 
is infected, each of the other healthy persons at that site 
has a probability r to become infected and with a prob- 
ability 1 — r it remains healthy. An infected person at 
time t becomes healthy at the next time step. For a cer- 
tain average density p, the average fraction of infected 
persons in the system fluctuates but maintains a steady 
time independent average A(r,p). In Fig. 4 we show 
that the average activity A(r, p) vanishes for r < r c and 
it continuously increases beyond r c . The threshold r c is 
the critical point of a phase transition from a completely 
healthy society to an infected society. The A(r, p) plays 
the role of the order parameter in this transition. We also 
notice that r c is in general a function of the population 
density p. In the inset of Fig. 4 we plot the variation 
of r c (p) with p. The value of the critical infection prob- 
ability decreases with increasing the population density 

1. e., more the density it is more likely that the infection 
really spreads. On the other hand, below a certain den- 
sity p < p c infection does not spread at all even with the 
maximum possible infection probability r c = 1. For the 
square lattice we estimate p c ~ 0.75. 

A number of different aspects of this model may be 
of interest. On average a human being remains more 
or less localized up to his/her home, home city or home 
country. Therefore perhaps it would be better to consider 
their motion as sub-diffusive (R 2 (t) ~ t 2fl with p < 1/2) 
rather than normal diffusion. Secondly DLFN may be 
important to study the reaction kinetic networks of two 
species reversible or irreversible chemical reactions A + 
B <-> C. 

To summarize, we have considered the evaluation of 
the mutual friendship network in a dynamic model of a 
society. Each member of the society executes a diffu- 
sive motion. Members of the society represent nodes of 
the network and their mutual friendships are the links. 
The dynamical random graph associated with the net- 
work shows a novel dimension dependent phase transi- 
tion. With a certain death/birth probability the net- 
work displays the six degrees of separation effect. We 
also observe that such a society remains always healthy 
if the average population density is below certain thresh- 
old value, which should have very important practical 
consequences. 
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